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Abstract 

The aim of this note is to study the dynamo properties of the magnetohydrodynamics 
system at vanishing 7?„,. Improving the analysis in [61, we shall estabUsh a generic Lyapunov 
instability result. 



1 The induction equation 

The incompressible Magnetohydrodynamics (MHD) system is a classical model for conducting 
fluids. It is derived from both Navier-Stokes equation for the fluid movement under the magnetic 
field and the Maxwell equations with Ohm's law for the evolution of the magnetic field generated 
by the fluid movement. It is in dimensionless form: 

dfU + u - Vu + Vp - = {V A b) A b + f, 
dtb -V A iu Ab) - J-Ab = 0, 
V ■ u = 0, 
^ V b = 0, 

where Re and R„, are the hydrodynamic and magnetic Reynolds numbers. In this system, u = 
u{t,6) and b = b{t,9) are respectively the velocity field of the fluid and the magnetic field, 
vector- valued functions in M?, and / is an additional forcing term (anything that is not the elec- 
tromagnetic force). To lighten notations, we will assume in this note that Re = I and R,„ = s, 
where e is a small parameter. Note that the Reynolds magnetic number - that is the one that 
corresponds to the definition commonly used - is defined by R,,, = where L and V are re- 
spectively the typical variation length and amplitude of the velocity field u. This means that for 
a given time-and- space scale for the velocity field, we consider large magnetic diffusivity (or 
vanishing conductivity) of the fluid. The reader may refer to ^ to have a complete explanation 
on how the above MHD system is computed. One may say it is the base system for describing 
MHD in the case of a conducting viscous fluid. 
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Our aim here, related to dynamo theory is to deal with the stability of solutions of this last 
system, and more precisely with the instability of solutions of the form 

{u,b) = {u{t,x),0). (2) 

The physical problem behind the study of the stability of such solutions is to understand the 
generation of magnetic field by the fluid. The fluid "gives" energy to the magnetic field by the 
term V a {u a b), which should compensate the dissipation term in the equation. We know for 
long that this kind of phenomena occurs in planets and stars like the Sun or the Earth and make 
them act like a magnet (See ^ where the problem was first posed). 

There have been many results in the past century in dynamo theory. First results where 
"negative" results in the sense that they gave condition for which fluids cannot generate a dynamo 
effect, i.e. this corresponds to situations where solutions of the form © were actually stable. 
Those results spread the idea that the velocity field has to be chaotic enough to expect a dynamo 
effect (See M.M.Vishik, [fT6ll ). The knowledge in the domain concerning stability of the MHD 
system has been summed up in |I3. 

Then positive results began to emerge concerning kinetic dynamos, that is a dynamo where 
we consider only the induction equation in the system: physically, this corresponds to the case 
where the Laplace force is neglected in the fluid movement, i.e. the magnetic field has no retro- 
action on the fluid. This case is much simpler since it corresponds to only a linear case. More 
recent results gave positive result with the full MHD system, which is non-linear. 

Li this note, we will focus on a particular mechanism for the generation of magnetic field, 
that is the "alpha-effect", which was first introduced by E.N. Parker [[TT|. This mechanism is 
based on scale separation: we formally decompose the fields into two parts, a fluctuating one, 
evolving on small lengths, and a mean one evolving on much bigger ones. The idea is that the 
small-scale field can have an effect on the big-scale one through the mean part of the crossed- 
term V a (a a Z?) in the induction equation and may create instability on the large-scale. Note that 
this kind of mechanism has been experimentally confirmed much recently by R. Stieglitz and U. 
Miiller lfT4l . 

Three regimes of R,„ can be identified in the MHD system: 

The intermediate regime where R,„ is of order 1. We considered this regime in a former arti- 
cle 0. 

The regime of large Rm. In that case, one can expect a dynamo at the scale of the velocity field, 
since the diffusive effect is very weak at such scale. A related question is whether the 
dynamo is "fast" or "slow", that is whether the best growing mode decreases to as R,„ — > 
+00 or not. So far this regime has been mainly studied formally and numerically (up to 
Rm ^ 10^) in physics lfT0l[T3l[T2l [T]|. in particular in the case of ABC-like flows, and the 
result seems to be difficult to predict in the considered cases. This regime will be the object 
of a forthcoming paper. 

The regime of small 7?,„. This regime is relevant to laboratory experiments, and is the focus 
of the present note. More precisely, we use the a-effect mechanism to obtain a linearized 
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instability result for system©, around steady states {u,b) = {U{6),0). Then, we conclude 
to a Lyapunov instability result for the full non-linear system. 

The base flows we consider satisfy 

Ue'H'^iTy (3) 

and 

f/ = 0, and Vg-U = 0. (4) 



We denote by ^ the set of fields U of this form. ^ is a Frechet space for the associated semi- 
norms 

\mi = j]\^num'- (5) 

We have 

u = U + u', b = + b', (6) 

focusing on the growth of the perturbations u', b' . As we shall see, it is also convenient to do a 
rescaling in time: t' = e^t. Dropping the primes, we end up with: 

' dtb - s^V a{U Ab)- s^V a{u Ab)- E-^Sb = 0, 
s^dtU + u ■ VU + U -Vu + u -Vu + Vp - Au = {V A b) A b 
Vb = 
V ■u = 

As mentioned before, a key point is the understanding of the linearized induction equation, 

dtb - s^V a{U Ab)- e-^Ab = 0, (8) 

Motivated by the analysis of the a-efi'ect, we will look for an exponentially growing solution of 
that equation based on scale separation: 

b%t, e) = e^^\b{s^e) + sb{s^e, e)), (9) 

where b = b{x) and b = b{x, 9) depend on the original small-scale variable 6, but also on a large- 
scale variable x = s^d. In this decomposition, b e 'H^^CM? x T^) is periodic and average-free 
with respect to the second variable 6, whereas b s ^-{^^^{M?). 

Using the formal variable separation above, the equation then writes under the form 

sA'^b + A'b - v., a{U Ab)- -V,,- a{U aI)- a{U Ab) 

1 1 ^ ^2 ' ^^^^ 

- -V, A (t/ A ^) - -Ash sA^b A,b - - 2 d,,e~b = 0. 

i=l 

We split this equation between mean part and fluctuating part: 

Agb + sVg A [U Ab) +Vg A {U Ab) + S^{V^ -Ve + Vg- V.0^ + 

e^V., a{U Ab)+ £^V,- A + s'^A.b = e^A'b, 



(11) 



v., A f/ A S + ^J) = X^b. (12) 
And we rewrite the fluctuating part of the equation to separate terms with b and b: 

Ae5 + sVe a ([/ a + s^{V, • Ve + Ve • V,)b + s^V , a (t^T^^ + e^^,b - e^A^b ^^^^ 
= -Ve a{U Ab)- e^V,- a {U a b) 

Equations (fTTj) to (fT3l) will be used extensively throughout the note. 

We can now give the theorems that we will prove in this note. The first one is related to the 
induction equation ([8]): 

Theorem 1 Let s > be a real number. There exists a dense subset Q. of 0^ such that for all 
U s Q., there exists a real Sq > 0, a vector ^ e and a function Au : [0, eo['^ R+ analytic such 
that for all s e [0, s^, there exists a solution b of ^ of the form: 

b\e, t) = e^"^'^'+'''^-\b' + sb'{e)), (14) 

where U e ^,b' eR^ and b' e 'H^{T^). 

Our second theorem concerns the full MHD system ([7]). For this next theorem, we introduce for 
any vector T e (R* )3 the Sobolev space "K^ = "H^R/TiZ x R/TjZ x R/T^Z) 

Theorem 2 Suppose s > 3/2 + 1. Then there exists a dense subset Q. of ^ such that for all 
U E Q, there exists a vector T = {Ti,T2,T^) with integer coordinates and a sequence {sn)n>o 

decreasing to such that the solution \ ^] of the (MHD) system ([T]) is nonlinearly unstable in 



IHj in the following sense: 
There exists initial values 



e -H^, (15) 



and Co > such that for all 5 > Q, the solution [^,^\ of ^ with initial value 5 \ satisfies 



' its) 



b. 



bsj \bo, 

> Co (16) 



for some time t^. 



This result extends substantially the article ^ by D. Gerard- Varet. Recast in our variables, 
this article considers the case of large magnetic diffusion, that is = e « 1 . Let us point 
out that D. Gerard- Varet didn't use the physicists scaling for the velocity field, thus making 
possible confusion to the fact that he chose = 1" and that it actually corresponds to the case 
where Rm is vanishing with s. To get the usual Reynolds number you have to do a time-space 
transformation t i-^ t/s^, x t-^ x/s^ (and multiply the amplitude V by 1/e), thus making the 
Reynolds number s instead of 1 . 
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Briefly, it is shown in (611 that for all m e N, there are solutions of ([8]) (on a larger box 
depending on e) that go from amplitude e'" initially to an amplitude 77 = rj{m) > independent 
of s (but depending on m). Note that, as j] depends on m, the result in does not yield Lyapunov 
instability at fixed small e, contrary to ours. 

The main idea to prove our two theorems is to construct an exact solution (b, b) of equations 
(fTTI) to (fT3l) . We shall obtain it through a perturbative argument, starting from the case e = 0. 

2 Case £ = 

2.1 An eigenvalue problem 

The fluctuating part in the case e = is 

= -V^, a{U A b). (17) 
Since the Laplacian is invertible as a function from "KJ^^ to IH^, we can write 

b = -A-'Ve a{U Ab) = ^{e)b, (18) 

where ^{d) '■ 'Hq '^q^^ can be taken as a parameter since it is not involved in 

Now we can replace b by ^{Q)b in the mean part of the equation (fT2)) to have an equation 
involving only b: 

v., A U A ^{e)b + A,b = A%. (19) 
We introduce the matrix a defined by 

a{U)b = U A ^{6)b = - | i7 a A^' Vg a ([/ a ^) dO. (20) 

(The name of the matrix is the one related to the so called alpha eff"ect mentionned in the intro- 
duction) 

Now the equation is 

v., A {a{U)b) + A„b = A%. (21) 

We will find a solution of this last equation under the form b{x) = e'^'^b'. Using this stanza, 
the equation then has the form (dropping the prime to lighten notations in the end of the subsec- 
tion): 

A {a{U)b) - \^\^b = A%. (22) 

Thus we have come now to an eigenvalue problem with the matrix ^ a a{U). The properties 
of the matrix a are the source of the alpha-effect mechanism; they have been already studied in 
Il6l and in ^ with a more general case. We give here a proposition which is only the particular 
case of the similar study done in those papers: 

Proposition 1 The matrix a{U) defined above is real and symmetric. Furthermore, there exists 
a dense subset Q. of such that for all profile U e Q, there exists some ^ e such that the 
matrix defined by the left hand side of (l22l) admits an eigenvalue of positive real part. 
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The reader should refer to those paper for a proof of this proposition. 

From now on, we will assume that our base flow belongs to this subset Q.. This gives the 
existence of a vector ^ g and an eigenvector bo associated to an eigenvalue /l° with positive 
real part such that 

b{x) = e'^'bo (23) 



is a growing mode of (121)) 

3 Existence of growing mode for £ > 

In this section we will search for a solution of the induction equation for e > 0, and we will 
search a solution under the form 

b%e, t) = e'^^''+'''^ %b + sb{e)), (24) 

based on the result of last section, where ^ is the one chosen for e = 0. Thus the equations 
become: 

^eb + eVe a [u aI) + is\^ -Vg + Vo- ^)b + eH^ a - e^\^fb - e^A'b ^^5) 

= ^Vgi A {u Ab) — sH^ A [u A b). 



AUAb- \^\^b = A'b. (26) 

The left part of (|25l ) may be seen as an operator from 'HJ'*'^ to "Hj, and for small enough s 
and A'^ close enough to /l" it is a perturbation of the Laplacian. Thus this operator is invertible 
and we can write again 

b = ^^,,AO)b, (27) 
where ^^^e^a^{9) = — L~]^^ o (Vg au a ■ + eH^ a (m a •)) is the operator from 'Hg to "TYq^' where 

L^,s.A^ = + eVe A (m A •) + • + • ^) + sH^ a (^7^) - - s'^A' (28) 
and the mean-part equation is 

A ia^^,Ai^)b) - \^\^b = A% (29) 
where the matrix a depends here on some parameters: 



a^,E,A^{u)b = u A ^^^s,A'^{9)b 

= U A Ljj^^^{Vg A U A b + sH^ A [u A b)) 

Jmv 



(30) 



/[0,1]3 

This equation admits non-zero solution b if and only if 

det(i^ A a^,sAu) - I^Pld -A^A) = 0, (31) 
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where 

A ag^^j.{u) = i \ ^3 -^i \a£^s,A'^{u) = A^a^^sj'^iu). (32) 

V-^2 fl 

f being fixed by the former section, we define the following function 

f{e,iu) = det(/f A a^,,,^(i<) - \^\Hd-fxld), (33) 
defined in a neighborhood V of (0, A^) in R^. 

4 Implicit function theorem 

We know from former section that /(O, /l") = 0. Now we want to prove that for e > small 
enough, there exists A^ such that f{s, A^) = 0. Furthermore, / is analytic in (e,//) in V. 

0, so that we can apply the analytic implicit functions 

theorem and conclude that there exists a neighborhood V7 of and a function ip: W ^ R such 
that ^(0) = /l° and Ve e W, f{e,tp{s)) = (we can then restrict this neighborhood to positive 
values of s). 

We compute now 1^ : using usual computation rules on composition of functions, we 
have that 



We prove in this section that ^ 



df 



e=0 



ddet|,.^^„^^^^(„)_l^pid_^oid o d(if A a^,s.M - 1^1 ^ -^lU)\e^o^ (0, 1). (34) 



The computation of the second differential is straightforward and gives 



d(if A a^,,,^(i<) - |f|^Id-//Id)|,=o (0, 1 



= - Id. 

£=0 

fj=A« 

(35) 

And the computation of the first one gives 



ddet|,^,,^^^„(„)„l^pid_,oid (-Id) = -Tr(com(if a a^^oA>^) - |f ^ Id Id)). (36) 

We know from the second section the matrix A = if a a^ ojo{u) — |f p Id ~A^ Id is diagonalizable 
and that his eigenvalues are 0, — /l", —2/1°. Thus, com(A) is similar to the matrix 

^0 

com I -i" 1=1 1, (37) 



.0 -2A 



which is of trace 2{A^Y ^ 0. 
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Thus 



5/i 



e=0 



2(^)2 ^ 0. 



(38) 



Since (/? is continuous on 0, and i^(0) = has positive real value, this remains true for s > 
sufficiently small, which concludes the proof of the existence of a growing mode for a small 
e > 0. 



5 Nonlinear instability 

We will now prove nonlinear instability of the system ([8]). For that we will adapt a method 
developped by S. Friedlander, W. Strauss and M. Vishik in (Si (itself adapted from Y. Guo and 
W. Strauss in jHl) to prove nonlinear instability given linear instability. The condition s > 3/2+1 
in the main theorem is used only in this section to obtain an inequality on the nonlinear terms 
similar to that in the article [5]. Actually, using the parabolic feature of the MHD system, it 
would be possible to lower our regularity requirements (see for instance Friedlander et al [H 
in the context of the Navier-Stokes equations). We stick here to regular data for simplicity of 
exposure. 

Let IJ ^ 3^ belonging to the dense subset of linearly unstable flows described in the previous 
section. 

We introduce Leray operator P to get rid of p. We rewrite system ([7]) as 



a,'^ +L,'^ = 2('^,'^), (39) 



where 



t.i^ = .-f''";"'t«'J"):,'^"') (41) 

' \ -V A {u Ab) - J ^ ^ 

is the linearized operator around , and 

is the nonlinearity. 

We shall consider this system in l-lj, T to be specified later. We also denote = Hj the 
space on the torus. Since we will work in these spaces until the end of the proof, we will drop 
the T in the notation. We have the following a priori estimates for the nonlinear term: 

Every term of every component of can be written as a product of one component 

of ^ and one of (composed by P). Thus, we have the first estimate 

||2(^,^)||l2 <.Cs-^\\\^\\V^\y ^Ce'^\\^y\\V^\\L^. (43) 
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For r = j{s + n/2 + 1) > n/2 + 1 (here n = 3), we have that 

II 11^30 < 11^ 11^,-, (44) 

and 

II II < C II ^11 ^,11^11^7 (45) 

for T] such that r = (1 - 77)5 (i.e. V = ^ ~ = ^ - j;)- 
Thus, we have: 

112(^,^)11^. < Ce-^||^||^+''||^||;7,'' (46) 
We shall prove the following theorem, which is a mere reformulation of the second theorem: 



Theorems Suppose s > 3/2 + 1. Then there exists a vector T = (Ti,T2,Tt,) with integer 
coordinates and a sequence {sn)n>o decreasing to such that the system (l39l) is nonlinearly 
unstable in 'Hj in the following sense: 

There exists a growing mode = ( f ° ) with ||'^o||'H' = 1 of the linearized system and a 

\boJ 

constant Co > depending only on '^q, s,n,p such that for all 6 > 0, the solution of the full 
magnetohydrodynamics system with initial value 6^q satisfies \\^{tg) ^ C^for some time t^. 

Remark 1 This theorem only proves that the solution reaches high values infinite time no matter 
how small the initial condition is. In particular it doesn't say wether it explodes infinite time or 
is defined for all time, since both situations are possible in the result of the theorem. Actually, for 
small enough initial data ( that is small enough 6 ), it is classical that the solution of the system 
doesn't explode infinite time. More details in / [731/ (Theorem 3.8) 

Assume in contrary that the system is nonlinearly stable while spectrally unstable, that is that 

for any growing mode of initial value = ( ) with ||'^o||'H^ = 1 in any box T of the system 

\boJ 

(i.e. an eigenvector of associated to an eigenvalue with positive real part) and for any £ > 0, 
there exists a S > such that for all t > 0, the solution (t) of the system with initial values 
6% satisfies ||'^(OI|-w^ < e. 

5.1 Proof of the nonlinear instability 

Let bii^O) = e"^'^'^b'^{6) be the unstable eigenmode given by Theorem 2. We fix e < sq in such a 
way that 

/ 2n 2n 2n \ 

belongs to N^. This is possible, as it is clear from the proof of Theorem [2] that ^ can be chosen in 
(2;7rQ* )^ (take s to be the inverse of a large integer, which can be chosen as big as we want). In 
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particular, '^o = ^ is an eigenvector of associated to an eigenvalue with positive real part 
in-K^. 

Denote by p the maximum real part of the spectrum of L^: 

p = max {Re A, As Spectr(L,)} . (48) 

By the previous remark on we know that p > 0. Moreover, since the spectrum of is 

only made of eigenvalues, there exists an eigenvector = f l° ) with eigenvalue A satisfying 

\boJ 

exactly Re /I = p. 

Moreover, by standard properties of the spectral radius of we know that 



(f = lim lle'^i'/'. (49) 

Thus, for all r] > 0, there exists a Cr, such that for all f > 0, 

lle'^'ll < C,e^('+?> (50) 

Then, the solution of (|39l ) with initial data S'^q can be written 

^ = ^ + (Je-^'^o, (51) 

where ^ satisfies 



= L,^ + !2(^,^) 
^f=o = 

Using Duhamel's formula, we can write the solution under the form: 



(52) 



^{t) = ( e^'-'^^^Q{'W,'W)ds. (53) 
Jo 



Define 

1 



= supjr > 0,\fs < t,\\^{s)\\L2 < -6e^'\\%\\^2^, (54) 

where A is the coefficient in the growing linear mode. Then, since we supposed that 'W = ^/ + ^ 
remains bounded, and \\%\\ ^ 6eP'\\'Wo\\, we have that Ts < oo, and for all t < T^, 

\\'^{t)U2<^6e'"\\%U2. (55) 

Using a priori estimates stated earlier, we have for 77 = ^ — ^±2 = 5 ~ 

< C,s-' ( eP^'+'^^^'-'^\\^\\[t'\\'^\\lj^.' ds. (56) 
Jo 
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Jo 

< Cr,6'^'>e'^'+^>s-^-''—{e'^-P - 1) ^^^^ 
VP 



Thus, getting back to : 

W'^y > H%\\L^eP' - C'(5^+''£-2-V('+'')^ (58) 
By definition of Tg, we have: 



(59) 



that is 

C'6'^s-^-\'^''^p, (60) 
^ in " " , — . (61) 

Denote by ts the quantity 

^ (2cJ-°'ia + ^)) ^^^^ 

Thus ?5 < Tg, and at this time ts 

„^,U....-^(„^.„.-^)^^„..,U.(M^)'". (63) 

Thus 11'^ (and \\'^\\^^) can be bounded below by a value independent of 6, which contradicts 
the assumption, and proves the theorem. 
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